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ON INTEGRATION OF HAMILTON-JACOBI 
PARTIAL DIFFERENTIAL EQUATION 


Abolghassem Ghassari 


SUMMARY 



The Hamilton-Jacobi partial differential equation is fundamental in plan- 
etary and lunar theories. The solution of many perturbations theories and also 
numerous problems of the control theory is reduced to the problem of solving 
the Hamilton-Jacobi equation. 


The application of the Hamilton-Jacobi theory to dynamical systems is based 
on the assumption that the canonical equations can be solved by Jacobi's theorem, 
which requires essentially the knowledge of any complete integral of Hamilton- 
Jacobi equation. 


The purpose of this paper is to illustrate a method leading to a complete 
integral of Hamilton-Jacobi equation, which is based on the variational principle 
and transversality condition. An application of this method to a particular case 
of Hamilton-Jacobi equation is described, and the classical example of wave 
fronts or parallel surfaces is considered. 
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ON INTEGRATION OF HAMILTON-JACOBI 
PARTIAL DIFFERENTIAL EQUATION 


INTRODUCTION 

The equations of motion of a system of n mass -points in terms of general- 
ized coordinates are given [l] 1 by Lagrange's equations. These equations are 
n differential equations of the second order with n unknown functions. 

The Lagrangian equations of motion have not been developed with the special 
reference to the problems of celestial mechanics, and are not yet adapted to the 
needs of the space mechanics. They have to be transformed into a system of 2n 
first-order differential equations called Hamiltonian or canonical equations of 
motion. 

In canonical form the 2n variables appear separated into two sets, essen- 
tially different in character, namely, into n coordinates of position and n co- 
ordinates of momentum. The canonical equations of motion are of the first order 
and therefore more tractable than the Lagrangian equations. These equations do 
not facilitate the solution of particular problems and there is no technique known 
for solving them in general; but, they form the basis for most theoretical dis- 
cussion since they are amendable to the very considerable body of theory that 
has been built around first order equations. 

The reduction of the order of a differential canonical system can, in theory, 
be performed by obtaining first integrals of the system. 

There are many classical methods such as Hamilton -Jacobi theory, Lindstedt's 
method, Whittaker's method, Delauny's method, the von Zeipel's method, and others 
for reducing the order of canonical systems. 

The most general technique which has been applied quite often to the prob- 
lems of celestial and quantum mechanics is the Hamilton -Jacobi method which 
consists in obtaining a contact-transformation 2 such that the new Hamiltonian is 
identically zero, and therefore the new variables are all constants. 


figures in brackets indicate the literature references at the end of this paper. 

A 

^For a rigorous and general principles of classical dynamics refer to Whittaker’s Analytical Dynam- 
ics [ 1 ] Chapters X-XII. 
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The integration of the canonical equations of motion can be made to depend 
on the solution of a first-order partial differential equation called Hamilton- 
Jacobi (H-J) equation, which does not contain the unknown function (generating 
function) explicitly. 

The integration of canonical equations is thus reduced to the integration of 
a first-order partial differential equation, which in no way means a simplifica- 
tion of the problem. For the treatment of partial differential equations is usually 
more complicated than that of ordinary differential equations. But, due to the 
important applications of H-J equation in dynamical problems and also in the 
theory of optimal control, the discussion of H-J equation has been very success- 
ful in illustrating the intrinsic features of the problem concerned. 

The canonical 2n equations can be integrated by Jacobi's theorem which is 
based on the previous knowledge of any complete integral of H-J equation, a solu- 
tion containing n arbitrary constants in addition to the additive constant. The 
integration of canonical equations is, therefore, connected to the general theory 
of partial differential equations of the first order. One realizes at once that when 
Hamiltonian function is time- independent, the energy integral constitutes a 
first integral, and the canonical equations are the differential equations of the 
characteristics of the 1I-J equation. 

General application of H-J method to planetary and lunar theories is based 
on the assumption that the canonical equations can be solved by Jacobi's theorem, 
which requires essentially a complete integral of H-J equation. To obtain, on the 
other hand, a complete integral of H-J equation one has to set up, formally, its 
characteristic system of differential equations which is the same as the canonical 
system of equations. Thus, it seems there is a plausible paradox. Nevertheless, 
Jacobi's method has been very helpful in many dynamical problems, for one can 
find out quite often, by some direct procedures and separation of variables, a 
complete integral of H-J equation without using the canonical equations. 

Hamilton-Jacobi equation can be integrated, in general, by the application of 
Charpit's method which provides an infinite number of complete integrals. Never- 
theless, the method of the separation of variables, subject to Levi-Civita's condi- 
tions, can also be applied. 

The purpose of this paper is to illustrate a different method leading to a com- 
plete integral of H-J equation, which is based on the variational principle and 
transversality condition. An application of this method to a particular case of H-J 
equation is described, and the example of parallel surfaces or wave fronts is 
considered. The application of this method to some problems of space mechanics 
and optimal control theoi’y will be treated later on. 
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METHOD OF INTEGRATION 

Extermals with Variable End-Points on Fixed Curves 
Let us consider the functional 


i[e ab ] = f f ( x ’ y. z ’ V> *> dx (!) 

J A 

taken along a variable admissible arc E AB (Figure 1) whose end-points A (a, c) 
and B (b, d) describe the two fixed and given curves C and D. 

Dot denotes the derivatives with respect to x. The coordinates x, y, and z 
along the arc E AB may be taken as functions of the parameter time t. 



Figure 1 


The necessary conditions [2] 3 that the admissible arc E AB extremizes the 
functional I[E AB ], i.e. 


SI EE ab ] sO, 

where S indloates the first variation are that: 
*Fer mere details, see Bill* [2] , Chapter VI. 
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1. The admissible arc E.„ represents a non-singular extremal of the func- 
tional I, i.e., y(x) and z(x) must have continuous first and second order deriva- 
tives on (a,b), and E AB satisfies the two Euler-Lagrange differential equations: 


_d_ 

dx 


Bf 

3y 


11 




= o 


> 


( 2 ) 


_d_ /_3f 
dx \ 3z 


_ 11,0 




where the function f is supposed to be twice continuously differentiable with re- 
spect to each arguments x, y, z, y and z. 


2. The admissible arc E AB satisfies the transversality condition 


or 



a a « n 

dx + dy + — dz = 0 

By Bz 


(3) 


f , 


dy\ Bf 
dx ) By 


(z -ll) 11 

\ dx / Bz 


at its intersections point A and B with the fixed curves C and D, i.e., the non- 
singular extremal arc E AB must cut the curves C and D transversally at the end- 
point A and B respectively. 

Needless to say that the arguments x, y, z, y, z belong to the extremal E AB 
and dx, dy and dz represent a direction tangent to the curves C and D at the points 
A and B respectively. 

The general solution of Euler's equations (2) form a four -parameter family 
of extremals. In order to determine these four parameters one has to verify that: 

1. The general solution satisfies the two unknown abscissas a, b of the vary- 
ing end-points A and B. 

2. The transversality condition is satisfied at both end-points A and B. 
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Elimination of a and b between the six conditions obtained leads to four rela- 
tions for the determination of the four parameters. 




The Euler equations (2) are of the second order and are, therefore, equiva- 
lent to the following system of the first order: 

dy . dz ''l 

— — = y > — = z 

dx dx 


B 2 f 

B 2 f . 

+ — ; — y + 

B 2 f . B 2 f 

z + 

dy ^ B 2 f 

dz 

=- ^ 

By Bx 

By By 

ByBz By By 

dx ByBz 

dx 

3y 

B2f 

B 2 f . 

+ y 

B 2 f . B 2 f 

+ — ; Z + — r 

dy B 2 f 

di 

Bf 

Bz Bx 

Bz By 

BzBz Bz3y 

dx + Bz3z 

dx 



where y and z are considered as new unknown functions. 

In system (4) the arguments y, z, y, z figure in a nonsymmetric manner and 
the equations are not solved for the derivative of the all unknown functions. 

The application of Legendre's transformation [3] 4 


P 


Bf 

3y 







( 5 ) 


h(x, y, z, p, q) 


. Bf . Bf 

y — + z — r 

By Bz 


- f 




4 See [3] Chapter IV, and [4] Chapter II. 
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leads to a more satisfactory first-order system: 


dy _ Bh dz _ Bh 
dx Bp ’ dx Bq 


> 


dp 3h dq _ 3h 

dx 3y dx 3z 


( 6 ) 


called the canonical differential equations of the variational problem, provided 
the Hessian of f ^ 0, i.e. 


D = 


3 2 f 
3y 3z 


t 0 . 


(7) 


Taking into account of Legendre's transformation (5) the differential of the func- 
tional I, taken along a variable non-singular extremal arc E AB whose end-points 
A and B describe two fixed curves C and D, can be written [2] 


dl = [-h(x, y, z, p, q) dx + p dy + q dz]® . (8) 

If we limit ourselves to the case where the extremal E AB is transversal only at 
curve C at A, then 


dl = [-h dx + p dy + q dz] . 

u 


(9) 


Suppose now by each point B of curve D passes one and only one extremal, 
then the value of I as well as the values of p and q depend on the coordinates 
x, y, z of variable end-point B. 

Equation (9) can be written: 

31 U, N 31 31 

— = - h(x, y, z, p, q) , = P- — = q 

3x 3y dz 
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or 


31 

3x 


+ h 



y. 




( 10 ) 


Therefore the functional I satisfies a first-order partial differential equation 
solved for 3 1/3 x, which does not contain the unknown functional I explicitly. 
This equation is called Hamilton-Jacobi (H-J) partial differential equation. 


Conversely, this result leads to a general method of integration of first- 
order partial differential equations not containing the unknown function. In fact, 
in solving such an equation for one of its derivative, say 3W/3x, it takes the 
form: 


3W 

3x 


+ H 



3W 3W \ 
z. — , — ) 
3y 3 z ) 


= 0 . 


( 11 ) 


The corresponding canonical equations are 


dy _ 3H dz _ 3H 
dx 3p dx 3q 


> 


dp _ £ dq = _ 3H 

dx 3y dx 3* 


( 12 ) 


where H = H(x, y, z, p, q). 

The integration of (12) gives y, z, p and q in terms of x and four arbitrary 
constants. 


Now we consider the functional 


I = 



f d x 
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I 


where f is the function of x, y, z, p and q defined by Legendre's transformation 
(5) and the canonical equations (6). Therefore 


3H 

f = p — + q 

3p 


3H 

Bq 


- H . 


(13) 


Substituting in f for y, z, p and q their expressions (solution of 6) in terms of 
x, we will get the functional I in terms of the coordinates of B, and it satisfies 
the H-J equation (11). Thus we obtain, by this method, a solution of (6) which 
depends on the arbitrary constants or functions, the coordinates of A or the 
functions determining the curve C. It is supposed, of course, that the equation 
to be solved may be considered as a consequence of a variational problem, that 
is to say, in applying Legendre's transformation to Euler's equation one can get 
canonical equations and vice versa. Therefore the resolution of the system 


dy _ 3H dz 3H 
dx 3p dx 3q 

for p and q is possible if the Hessian of H ^ 0, i.e. 


(11 L\ /JUL \ _ / _£h_\ ! i 0. < 14 > 

\3p3p/ \3q3q ) \3p3q/ 

This implies that the function H (Hamiltonian function) must be different from an 
homogeneous function of the first degree in p and q, which excludes the case f = 0 . 

As a result we can deduct that this method provides us the integrals of the 
given equation. In order to obtain a general solution of H-j equation one can take 
an extremal transversal at curve C at the end-point A. If the end-point A is taken 
as a fixed point, one obtains for the functional I an integral whose integrand de- 
pends on x and four arbitrary constants, the same number of arbitrary constants 
which arises from the integration of the canonical system. These arbitrary con- 
stants are the initial values of y, z, p and q for x = a. The initial values of x, y 
and z represent the coordinates of the end-point A. It remains two more arbi- 
trary constants which will be used for describing that the extremal passes also 
through the end-point B. 

Therefore, in varying the end-point A on curve C, one obtains an integral 
depending of these two arbitrary constants which are the coordinates of A and 
one additive constant y, since W + y is also a solution. 
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PARTICULAR CASE OF HAMILTON-JACOBI EQUATION 

The application of the preceding method to the following special case of 
Hamilton-Jacobi equation: 


^ tH (» 

3t \3x By/ 

gives easily the most general expression for its complete integral. 
In fact, setting 


( 15 ) 


we get 


3W 3W 

= p, — = q 

3x By 


f = P — + q — - H = R(p, q) (16) 

Bp Bq 

where R is a known function of p and q. 

The Hamiltonian function is only a function of generalized momenta p, q and 
the corresponding canonical equations have the two first integrals 

p = cont. □ a, q = cont. - b (17) 

The momenta are, therefore, preserved under a contact -transformation, and the 
cyclic coordinates x, y are expressible as linear functions of time: 


3H(a, b) . N a 3H(a, b) 
x - a - — L_1 — ' (t - t n ), y - /3 = — — : 

Ba Bb 


(t - O 


(18) 


where t 0 and t correspond to the fixed end-point A (a,/3 ) and variable end-point 
B (x, y) respectively. 
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The trajectories (18) represent the characteristics of (15) and also the ex- 
tremals of the functional 


I[E AB ] = f f dt = f R(a, b)dt = (t - t 0 )R(a, b) 

•'A 


If the Hessian of H ^ 0, i.e., 


#*.) (#± ) - (Jii-j * o, 

3a3a/ \3b3b I \3a3b / 


a and b can be obtained from the system (18). 


(19) 


In fact, setting £ = BH/3 a, v = 3H/3b we get, according to the properties 
of Legendre's transformation 


- df(^, rj) , _ Bf(£, rj) 


R(a, b) = f(£, rj) 


Therefore 


I = (t - t 0 )R(a, b) = (t - t Q ) f(i, v ) 


= (t - V f 


X - a y - ft 

t ~ t o t - t„ 


- F(t - t Q , 


which gives the most general expression of the complete 


x - a, y - (3 ) , 
integral in the form 


W = y +F(t - t 0 , x - a, y - p) 

where F is a certain homogeneous function of the first degree in t - t Q , x -a , 
y - (5, and y is an additive constant to W, since the equation still stands if we 
replace W by W +y . 
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The function F is obtained, through the known Hamiltonial H, by 


F(t - t 0 , x - a, y - fi) = (t - t Q ) 



and the relation (16). 

As the motion of the free particle governed by equation (15) is conservative 
its solution (complete integral) can be obtained by applying the method of sepa- 
ration of the variables. This method gives the complete integral 


W = -ht+ax+/3y+y (21) 

which is a special form of (20), where a, fi, y are arbitrary constants and h 
is another constant called energy integral. 


EXAMPLE 

Because of the analogy between the motion of particles and the propagation 
of light rays one can also apply the method discussed in this paper to the special 
form of eikonal equation. 


(V W) 2 = 



( 22 ) 


which is fundamental in geometrical optics. Each solution of (22) describes a 
definite beam of light rays, and the wave surfaces are the surfaces of constant 
eikonal, i.e. the parallel surfaces of the form 


W(x, y, z) = cont. 

Equation (22) has the complete integral 

W = ax + fi y + (1 - a 2 - fi 2 ) 1//2 z + y (23) 
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or more specifically 


W = [(x - a ) 2 + (y - fi) 2 + (z - y) 2 ] 1/2 , (24) 

which represents a family of parallel surfaces with a , /3, y and w as arbitrary 
constants. 

From the eikonal equation (22) is easily deduced that we are dealing with 
the ordinary Euclidean metric 

ds = (dx 2 + dy 2 + dz 2 ) 1/2 , 

the geodesics of which are of course the straight lines. 
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